q_hig_s=pars(1);
zlow_s=pars(2);
q_low_s=pars(3);
zhig_s=pars(4);
q_mid_s=pars(5);
zmid_s=pars(6);
rbar=pars(9);
r=pars(11);
delta_s=pars(13);
xi_s=pars(14);

coupon_s=pars(17);
rho_s=pars(18);
sigma_s=pars(19);
phi_s=pars(20);
a_s=pars(21);
varphi_s=pars(22);
nu=pars(23);
markupsub_s=pars(24);
m_s_ss=pars(26);
b_up=pars(27);
b_dn=pars(28);
b_md=pars(29);
B_s=pars(30);
theta=pars(32);
l_s=pars(33);

A=pars(34);
y=pars(25);
A_s=A;
y_s=A^(1/(1-nu));
z_0_s=pars(12);
k_s_ben=pars(7);
varzeta_s=pars(8);
tau_s_ben=pars(10);
chi_s=pars(16);
B0_s=pars(30);
R0_s=A/nu-chi_s;
zeta_s=pars(31);

%parameters that scale with value added
% A_cailb=pars(34);
% y_calib=pars(25);
% A_s=A;
% y_s=A^(1/(1-nu));
% z_0_calib=pars(12);
% z_0_s=z_0_calib*(y_s/y_calib)/(A_s/A_cailb);
% k_s_ben=pars(7)*(A/A_cailb)*(z_0_s/z_0_calib);
% varzeta_s=pars(8)*(A/A_cailb)*(z_0_s/z_0_calib);
% tau_s_ben=pars(10)*(A/A_cailb)*(z_0_s/z_0_calib);
% chi_s=pars(16)*(A_s/A_cailb);
% B0_s=pars(30)*(A/A_cailb)*(z_0_s/z_0_calib);
% R0_s=A/nu-chi_s;
% zeta_s=pars(31)*(A_cailb/A)^(1/theta);

